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We investigate an alternative to the Sequential Propagator Method used in Lattice QCD calcula- 
tions of semileptonic form factors. We replace the sequential propagator with a stochastic prop- 
agator so that, in principle, all momentum and sink smearing combinations are available with 
only a single spin-color inversion. Practically, the stochastic noise is significant and must be re- 
duced at the cost of more inversions. We study the behavior of the stochastic noise and compare 
the computational costs of this stochastic technique and the Sequential Propagator Method. We 
also present preliminary semileptonic form factor results using the stochastic technique on Nf = 2 
configurations with a non-perturbatively improved Sheikoleslami-Wohlert action generated by the 
QCDSF collaboration. At a fixed cost, measured in terms of the number of heavy-quark inver- 
sions, the method provides more correlators for the extraction of the form factors at various q 2, s 
than the Sequential Propagator Method. These additional correlators reduce the total statistical 
errors of certain kinematic points, although the stochastic error is still comparable to the gauge 
error at other points. 
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1. Introduction 

Experimental measurements of heavy-light semileptonic decays, combined with theoretical 
input, can be used to extract the Cabbibo-Kobayashi-Maskawa (CKM) matrix elements \V UD \, \V CD \, 
\V c d\, and \V CS \. The determination of these matrix elements provides constraints on the CKM 
Unitarity Triangle and thus test the Standard Model. Conversely, |V C( /| and |V C . S | are known with 
high precision, and can be used to test the Lattice techniques used to calculate the decay rates 
of D mesons. The approach to current Lattice calculations of the semileptonic decay rates of D 
mesons involves constructing the appropriate three-point function using the so-called sequential 
(or extended) propagator. We investigate an alternative method using stochastic techniques, which 
we refer to as the Stochastic Sink Method (SSM), in the hope of achieving an overall savings in 
computational effort. With this method all momentum and sink smearings are in principle available 
using only a single spin-color inversion. In this report we present a basic comparison at fixed cost 
of the two methods, and provide preliminary results of the form factors using the SSM. 

In the following we focus on the decays of a heavy-light pseudoscalar charm-like meson (H = 
D) to a light-light pseudoscalar meson (P = 7i/K) and leptons (Z,V/). For these processes the 
differential decay rate can be parametrized as 

* m ^ m =^J^*w)r +{ w, <>■» 



S form factor 



Perturbat ively known 



where q 2 = (pn — pp) 2 is the squared difference between the initial and final state four-momentum. 
The greatest source of uncertainty in the theoretical calculation is due to the non-perturbative in- 
teractions parametrized by the form factor F + (q ). 

These interactions appear in the hadronic matrix element (P(pp)\V^ (q 2 )\H(pH)}- = WcYn Wl 
is a weak flavour-changing vector current, where yfc i s the charm quark and yfl is the d or s quark. 
The matrix element can be parametrized as a linear combination of the form factors F + and Fo, 

(P(pp)\V^q 2 )\H(p H )} = {p H + p P -q(m 2 H -m 2 P )/q 2 }^F + (q 2 ) + {q(m 2 ^ 

(1.2) 

On the lattice the matrix elements are extracted from three-point correlators with the following 
form, 

C 3 (T,t;p H ,q) = Y^-^e^^YsWcixJ) ■ fc^Vi&O ■ i/W H (0,0)|0> = 

x,y 

-£^"V ? Tr(^^^^ (1.3) 
where \j/ u is the spectator light-quark and M~ l is the propagator for quark x. In the limit of large 



time separation Eq. ( L3[ ) has the form 



T hmC,(T,f,PH,q) - |^|| X (P(p P )\V^\H(p H )) x e^ 7 ^\ (1.4) 

so that a determination of the amplitudes and energies from ratios or simultaneous fits with meson 
propagators can isolate the matrix element. 
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2. Stochastic Sink Method (SSM) 

The standard method for calculationing D meson semileptonic three-point functions uses se- 
quential propagators. The sequential propagator provides a way to calculate a heavy-quark propa- 
gator that connects all spatial sites x at the sink time-slice T, to all sites y and t at the vector current. 
It starts by taking a single time-slice of the spectator quark propagator, M~ 1 (x , T ; , 0) . The desired 
sink momentum p H , sink smearing W s , and sink gamma Ft are then inserted to get the "sequential 
source". The heavy-quark action is then inverted on this "sequential source" to get the sequential 
propagator, 

XX? 1 (y,t;x, T)e ip »*r f W s M- l (x, T;0,0). (2.1) 

X 

The sequential propagator can then be combined with the daughter light-quark propagator, M/(y,f;0,0), 
and appropriate gamma matrices to get Eq. (13). This method requires a heavy-quark inversion 
for each distinct sink momentum pn, sink smearing Ws, and sink gamma Ff. The computational 
effort required for this procedure can become prohibitive if many sink momenta and/or smearings 
are needed, as would be required in the use of the Variational Method [Q, |5|] for studying excited 
state decays. 

An alternative method which may be more efficient is to replace the sequential propagator 
with an all-to-all propagator We construct all-to-all propagators by generating random vectors 
Tjj r ' (x,T), r = 1, . . . ,7V, at a particular timeslice T using complex Z2 noise, with the property 

^ 1 4 ] (*) nf r] (z) = s x Aj + *(i/V3v) (2.2) 

where i,j label spin and color. We then invert the charm-like quark's Dirac operator, Mij(x,y), on 
each source r\ ■ (x,T) to obtain the solutions y/j (y,t): 

M kj (z,T;y,t)yf ] (y,t) = 4 ] (z,T) -> W f{y,t)=^M^{y^zJ)nl\zJ). (2.3) 

Z,k 

The average over the product of the sources and solutions provides an estimate for the all-to-all 
heavy-quark propagator, 

^ I ^\yM [r \x, T) = MjS&tix, T) +l t MJ k l (y,t;z, T) ^%T)ij^%T) ~ SgSu), 

ly r v ' %\ ly r 

all-to-all v v ' 

errors ^(1/VN) 

(2.4) 

where the stochastic error decreases with N, the number of source/solution pairs used. 

A stochastic estimate of Eq. (L3) can be constructed by combining the spectator and daugh- 
ter point-to-all light quark propagators and the all-to-all heavy-quark sources and solutions in the 
following manner, 



- Lr Tr (L* e-^TtM- 1 (0, 0;x, T)T f r] W (x,T) ■ J^e? V &t)TMi (y,t;0, 0) / 

C 3 (T,t;p H ,q) + 0(\/VN), (2.5) 

where the appropriate propagator smearings must be applied and = F, ■ = 75 and T = Jn. 
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Par am 


16 3 x 32 


24 3 x 48 


n cfgs 


100 


221 


a 


0.089 fm 


0.076 fm 




929 MeV 


272 MeV 


Wln,valence 


929 MeV 


438 MeV 


m D 


2.44 GeV 


1.93 GeV 



Table 1: Parameter details for the 2 ensembles used in the cost comparison. 

If the stochastic error term were negligible compared to the gauge noise we would have all 
sink smearings and momentum available with a single heavy-quark inversion. For the parameters 
in our calculation the error term is not neglible and must be reduced by additional noise vectors 
and/or noise reduction methods. Both improvements involve a computational overhead which must 
be accounted for in any comparison between the SSM and Sequential Propagator Method. We tried 
all combinations of spatial even/odd, color, and spin partitioning [g|], and found spin partitioning on 
its own to be the most computationally efficient method of noise reduction. The number of vectors 
used will be discussed in the next section. 



3. Sequential Propagator Method versus SSM and Preliminary Results 

We perform a simple comparison of the two methods making basic assumptions about the data 
set desired for our form factor calculation. We compare the total statistical errors of correlators 
constructed from the two methods at fixed cost, where the cost is measured in the number of 
heavy-quark action inversions. All available rotationally equivalent correlators corresponding to 
each individual q 2 point are averaged over to improve the statistics. The sink is placed at the 
midpoint of the lattice, allowing us to fold the data along the time axis. 

Two QCDSF ensembles with 2 dynamical sea-quarks ^ were used in the comparison with 
the parameters shown in Tab. [jj The ensembles were generated using a Wilson Plaquette action for 
the gluons and a non-perturbatively improved Sheikoleslami-Wohlert action for the fermions. Both 
actions have errors starting at 0(a 2 ). The interpolating fields in the correlators were Wuppertal 
smeared using the smearing factor and number of iterations that optimize overlap with the 
light-light meson ground state. It should be noted that the number of smearing iterations increases 
the stochastic noise significantly faster than it increases the gauge noise, and thus using a more 
customized smearing (with less iterations) for the heavy-light state H would reduce the stochastic 
noise presented in this report. 

We assume four sink momenta are desired: pn = (0,0,0), (1,0,0), (1, 1,0), (1, 1, 1). These 
sink momenta create data in the physical region, q 2 > 0, for the meson masses used in our calcula- 
tion. We also assume a minimal smearing basis of two different sink smearings, W s . The cost re- 
quired to build this data set using the Sequential Method is thus \2{spin/ color) x 4(pn) x 2(W S ) = 
96. The cost for the SSM with spin partitioning is 4(spin) x N, so using N = 24 is of comparable 
cost to the Sequential Method. Note that all pu are generated at negligible cost with the SSM. 

Figs. l(a)| - 1(d) are representative of the range of behaviour of the percentage statistical errors 



in our data set using ,/V = 24. Figs. 1(a) and 1(b) show the q = q max kinematic point for the 16 x 32 



and 24 3 x 48 lattices respectively. At this kinematic point the initial and final state mesons have 
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zero spatial momentum. The total noise for the correlator on the 16 3 x 32 lattice shown in Fig. 1(a) 



is dominated by the gauge noise-in fact this is true for all other kinematic points on this lattice. The 
stochastic noise for the 24 3 x 48 lattice shown in Fig. |l(b)| is dominated by the gauge noise near the 
heavy-light meson (t = 24), but becomes comparable to the gauge noise at the light-light meson 
(f = 0). 

The SSM correlator shown in Fig. 1(c) is constructed from an average of all available rota- 
tionally equivalent correlators, 6 in total, whereas the Sequential Method generated correlator has 
only a single correlator available. The averaging causes the SSM correlator to actually have smaller 
errors than the Sequential Method generated correlator. This suggests that there is statistical gain to 
be had in so many additional sink momenta. Fig. 1(d) shows the noisiest SSM generated correlator, 
relative to the Sequential Method correlator, that we found in our data set. In this case both the SSM 
and Sequential Method correlators can be constructed from 6 rotationally equivalent correlators. 

We have examined all correlators with an appreciable signal and draw the conclusion that 
with the modest number of 24 stochastic estimates the stochastic error is at worst comparable in 
magnitude to the gauge noise, while at best the additional data available results in smaller total 
errors for the SSM. We are expanding this analysis and expect to have more quantitative results in 
the near future. 

Bare lattice results using the SSM for F + (q 2 ) and Fo(q 2 ) calculated from the 16 3 x 32 ensemble 
are presented in Fig. ^[ 600 configurations with 24 stochastic vectors were used with all rotationally 
equivalent correlators averaged to improve statistics. 



4. Additional Considerations and Outlook 

The effects of &{a) improvement have also been examined and do not change our conclusions. 
The matching coefficient Zy is known non-perturbatively for these lattices ^ and the coefficient 
of the improvement term cy is known to one-loop perturbatively [|7|]. 

We've also investigated using the "one-end trick"[]|, ^] with one stochastic vector, with and 
without spin partitioning. In our examination, comparing the three methods with no momentum 
averaging, we saw a noise reduction for certain q 2 's in the one-end generated correlators. For all 
q 2, s however, after averaging over the available correlators for the SSM and Sequential Method 
(and no averaging for the one-end correlators), the one-end correlators' errors were larger. Because 
generating the additional rotationally equivalent correlators or stochastic vectors for the one-end 
method would make it clearly more expensive than the other methods, we've concluded the one- 
end method is less efficient for the parameters we are working with. 

The SSM is potentially computationally more efficient than using sequential propagators. With 
the SSM all sink momentum and smearings can be generated with a fixed number of heavy-quark 
inversions. Whether this method saves computational effort over the Sequential Method depends on 
the statistical improvement that additional sink momenta and smearings provide for the extraction 
of the form factors, and a more comprehensive investigation is underway. At this stage it can 
be stated that at a fixed cost the errors at particular q 2, s can be reduced by using the SSM. This 
method should also significantly reduce the cost required to use the Variational Method in three- 
point calculations, where multiple sink smearings are required for each sink momentum. 
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Figure 1: Four correlators corresponding to different q 2 kinematic points are shown. The upper plot of 
each figure shows the sequential (red) and stochastic (blue) correlators scaled by the appropriate two-point 
functions, such that the vector matrix element is proportional to the resulting plateau. The lower plot of each 
figure shows the percentage statistical errors of the two methods. The correlators shown are chosen because 



they are representative of the behaviour of the statistical errors in our data set. In Fig. 1(a) and Fig. l(b' the 
temporal component of the vector current at q 2 = q 2 mn for the 16 3 x 32 and 24 3 x 48 ensemble are shown. In 



Fig. 1(c) the temporal component with \pp\ = 1 and \pg\ = 1 is shown, where the stochastic correlator has 



been constructed from the average of the six rotationally equivalent correlators. The Sequential Method has 
only one correlator available in this data set. Fig. 1 (d)| presents the spatial component of the matrix element, 
where \pp\ = 1 and \pn\ = 0. Both the sequential and stochastic correlators have been averaged over the six 
available rotationally equivalent correlators. 
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Figure 2: Preliminary results using the SSM on the 16 3 x 32 ensemble with N = 24- and 600 configurations. 
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